Abstract. PSR B0943+10 is known to show remark- 
ably stable drifting subpulses, which can be interpreted 
in terms of a circumferential motion of 20 sparks, each 
completing one circulation around the periphery of the 
polar cap in 37 pulsar periods. We use this observational 
constraint and argue that the vacuum gap model can ad- 
equately describe the observed drift patterns. Further we 
demonstrate that only the presence of strong non-dipolar 
surface magnetic field can favor such vacuum gap forma- 
tion. Subsequently, for the first time we are able to con- 
. strain the parameters of the surface magnetic field, and 
I model the expected magnetic structure on the polar cap 
C\| . of PSR B0943+10 considering the inverse Compton scat- 

^ ' tering photon dominated vacuum gap. 
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1. Introduction 



The fundamental problem of pulsar research is lack of 
consensus about the radio emission mechanism, except 
that the electron-positron pair production is an essen- 
tial ingredient, to explain the observed coherent radiation. 
Moreover, the observed phenomenon of drifting subpulses 
serves as an excellent diagnostic tool to investigate the 
mecha nism of pulsar radio emission. PSR B0943+10 with 



period \P ~ 1.09 s and period derivative P = 3.5 ■ 10 
s/s exhibits long and stable sequences of drifting sub- 
pulses (Deshpande & Rankin 1999, 2001, hereafter DR99 
and DROl, Asgekar & Deshpande 2001). This pulsar is 
observable only at frequencies below 800 MHz (high fre- 
quency radiation most probably misses the line-of-sight 
due to specific observing geometry) and drifting subpulses 
with similar features are observed at 35, 111 and 430 
MHz, implying broad-band nature of the phenomenon. 
This suggests that drifting subpulses in this pulsar are 
associated with a stable system of emission subbeams, ro- 
tating around the pulsar magnetic axis. The aliasing re- 
solved analysis of the harmonic-resolved fluctuation spec- 
tra related to drifting subpulses yielded in precise esti- 
mates of periodicities of phase and amplitude modula- 
tion, namely P3 = 1.87 P and P3 = 37 P, respectively 
(DROl). These two periodicities are harmonically related 
and are used to determine the number of subpulse beams 
N = P3/P3 — 20 contributing to the observed modu- 
lations. The subpulse "drift-bands" clearly visible in the 
observed drift pattern (Backer 1973) are separated verti- 
cally by P3 (in pulsar periods) and horizontally by P2 (in 
pulse longitude). The latter is estimated as being about 
11° and can be converted into the magnetic azimuth angle 
7] separating adjacent subbeams. This requires knowled 



Karastergiou et al.(2001) observed simultaneously sub- 
pulses in single pulses of PSR B0329+34 at 1.4 and 
2.7 GHz. Careful analysis of these observations suggests 
strongly that a single plasma column fed by a single entity 
at the bottom of the pulsar magnetosphere is responsi- 
ble for the subpulse observed across the pulsar spectrum, 
which can be interpreted naturally in terms of sparks op- 
erating on the polar cap surface. This picture is probably 
more general and can be applied to other pulsars, includ- 



ing those with drifting subpulses. Deshpande & Rankin 
(1999D noticed that the subpulse drift phenomenon ob- 
served in PSR B0943+10 is consistent with the model pro- 
posed by Ruderman & Sutherland (1975, hereafter RS75), 
in which the subbeams are associated with sparking dis- 
charges operating within the vacuum gap (VG) developed 
just above the polar cap. The observed subpulse drift 
rate in this model arises from the speed of circulation 
of sparks around the magnetic axis, which is determined 
by rate of the E x B drift of electron-positron plasma 
confined to the isolated spark filament within the VG. 
If the planes of surface magnetic field lines are more or 
less perpendicular to the polar cap boundary, then the 
peripherial sparks will slowly rotate around the magnetic 
axis. Deshpande fc Rankin (1999 ) demonstrated clearly 
that the non-corotating subbeams associated with drift- 
ing subpulses actually lag the stellar rotation. This finding 
strenghtens the VG plasma E x B drift hypothesis. How- 
ever, the formation of VG itself suffered criticism owing 
to the problem of binding energy of ions, not being high 
enough to be retained on the stellar surface (Hillenbrandt 
& Miiller 1976, Kossl et al 1988). Morover, the association 
of drifting subpulses with sparks was questioned based on 
the fact that sparks will move towards the magnetic pole 
at a speed much faster than the E x B drift velocity, thus 



gp quenching t he stable drifting patterns (Gheng fc Ruder 



of the Inclination angle n- a,nd the impart a,ng1e 13. which 
for this pulsar is estimated from the polarization mea- 
surements and from the viewing geometry that leads to 
missing radiation above ^ 800 MHz. Using a = 11°. 64 
and (3 = -4.31° (DROl), one gets 77 = 18°, which in turn 
gives N — 360°/?/ = 20 subbeams rotating at the beam 
periphery of PSR B0943-f-10. 



man, 1977 ; Filippenko fc Radhakrishnan, 19821) . In this 



paper we demonstrate that indeed VG can form in PSR 
B0943-f 10 if the surface magnetic field is very strong and 
non-dipolar in nature. For the first time we are able to 
constrain the parameters of the surface magnetic field as- 
suming that the VG is dominated by the inverse Compton 
scattering photons. We also demonstrate that the embar- 
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rassing problem with rapid spark motion towards the pole 
will not occur due to specific geometry of the surface mag- 
netic field in the VG, at least in PSR B0943+10. 



2. Constraining surface field in PSR B0943+10 

2.1. Vacuum Gap formation 

Gil & Mitra (2001, hereafter GMOl) argued that VG can 
form in neutron stars with very strong non-dipolar sur- 
face magnetic field Bg — b ■ Bd, where Bd ^ 2 ■ 10^^(P • 
P_i5)°-5 G is the inferred dipolar field (B^ 4 x lO^^ G 
in this case) and P-15 = For surface mag- 

netic fields with 6^1, the radius of curvature of field 
lines TZ should be relatively small, with normalized val- 
ues TZq — TZ/R « 100, where R = 10^ cm is the neu- 
tron star radius. GMOl considered strong magnetic fields 
B > 0.1 Bq, where Bq = 4.4 x 10^^ G is termed as 
the quantum critical field. In such strong magnetic fields, 
high energy photons with frequency to produce electron- 
positron pairs near the kinematic threshold with photon 
energy huj — 2mc^ / smO, where sin 6* = Iph/TZ, Iph is the 
mean photon path for pair formation. The typical photon 
energy is Uoj — 1.5 ■ h'y^c/TZ in the case of curvature radia- 
tion (CR) seed photons (e.g. RS75) and hu — 2jheB /mc 
in the case of resonant inverse C ompton scattering (ICS) 
seed photons ( Zhang et al., 199^ ); see GMOl for details of 
the Near Threshold Vacuum Gap (NTVG) formation. 

The condition for VG formation in neutron stars with 
J7 • B < can be written in the form Ti/P, > 1, 



(5 X 103)7eg "&-ip~°-36p-{l.'^ cm (GMOl; their eq. 12) 
and the corresponding relation for gap formation is. 



where = 6 x 10'^b'^ '^^{P ■ P-is)"'^^ K is the criti- 
cal temperature above which ^^Fe ions cannot be b ound 
([Abrahams fc Shapiro, l"99l| ; |Usov fc Mehose, 1995]) . The 

/ • n1/4 

surface temperature Ts ~ Ik-eAVN/airrpj , where 
rp = 5-o.5-[^Q4p-o.5 jg actual polar cap radius, 
■^GJ = '^i^pBs/ (eP) is the kinematic Goldreich- Julian flux 
through the polar cap surface, AV — (27r/cP)Psft.^ is the 
potential drop across the polar gap. The efficiency co- 
efficient k in AV describes departure from black body 
conditions on the polar cap surface; k — qrad/lheat < I7 
where Qrad = crT^, queat = -{^ ■ Bs/27rc) • AV and 
cr = 5.67 X 10~^erg cm.~^ K~^s~^ . One can show that at 
the neutron star surface below the sparking gap the value 
of k can be as low as 0.3 (Gil 2002, in preparation). Consid- 
ering the CR seed photons as a source of electron-positron 
pairs one obtains the height of the CR dominated NTVG 
in the form hcR = (3 x 103)7^°■296-0■43po.2lp-a2l 

(GMOl; their eq. 6)|5 Using it for VG condition TjTs > 1, 
we get the relation between b and TZq in the form 
b = 15.6 X 7e°-3i (for k = 0.5), (1) 
and 6 = 20 X TZ^'^^ for k — 1.0. In the case of the resonant- 
ICS dominated gap, the height of the NTVG is hjcs = 



& = 5.1 X 7^g•38 
b = 6.4 X 7^g■38 



^ Note the typographical error in GMOl eq. 6 for powers of 
P and P which is corrected here. 



(for k = 0.5), (2) 
(for k = 1.0). (3) 
We argue in section 4 that from our modelling CR- 
NTVG is not likely to form in PSR B0943-hl0 and that 
for ICS-NTVG the value of the coefficient k is slightly 
lower than unity. Hereafter we only consider ICS-NTVG 
for further calculations. 

2.2. Complexity parameter 

Following the original idea of RS75, Gil & Sendyk (2000; 
hereafter GSOO) argued that the polar cap is populated 
with Njyiax ~ sparks with characteristic dimension T> ~ 
h, separated from one another by a distance d ^ where 
a — rp/h is the complexity parameter (see also Fan et 
al. 2001). Keeping in mind the grazing viewing trajectory 
in PSR B0943-I-10, we consider an arrangement of = 
20 < Nmax sparks on a periphery of the polar cap in 
this pulsar. Since the exact separation d ~ /i of adjacent 
sparks is unknown, one can write d — h{l + S) where 
\S\ ^ 1. We will arbitrarily choose the value of = 0.25, 
which determines the parameter space marked by the grey 
area in Fig. 1 (note that if \S\ — then the parameter 
space reduces to single line midway between lines (4) and 
(5), which is both unrealistic and inconvenient for further 
analysis). We can thus write N -V + N -d^ 20/i(2 + 5)^ 
27r(rp — h)ova — 7.34-1-3.18 5. Using the ICS dominated 
gap height h — hjcs we obtain the relation between b and 
7^6 in the form 6°■57^g "'^"^ = 2 + 0.86 5, which gives two 
limiting relations of the form 

6=4.9x7^^■l4 (for (5 = +0.25), (4) 
6 = 3.2 X 7^^■l4 (for 5 = -0.25), (5) 
corresponding to lines (4) and (5) in Fig. 1, respectively. 
Any point {b, Rq) between these lines corresponds to the 
VG with 20 sparks of diameter V = h, separated from 
each other by d = h{l + 5) (where -0.25 < 5 < -1-0.25), 
occupying a peripheral ring-like area on the polar cap of 
PSR B09434-10. Although the maximum possible number 
of sparks on the polar cap of this pulsar is Nmax ~ ~ 50 
(GSOO, Fan et al. 2001), only the 20 outermost sparks 
can contribute to the observed radiation due to the graz- 
ing line-of-sight geometry. Moreover, lack of conditions for 
vacuum gap formation towards the central parts of the po- 
lar cap can result in the hollow-cone structure of emission 
beams (see next section; Fig. 3). 

2.3. The'ExB drift rate 

The isolated columns of spark plasma within VG perform 
a slow E X B drift in the direction perpendicular to the 
planes of magnetic field lines, which in turn should be per- 
pendicular to the polar cap edge. Following RS75, one can 
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Fig. 1. Plot of the normalized surface magnetic field h = 
Bs/Bd versus normalized radius of curvature TZe = 71/ R. 
Lines (1) to (6) in the above figure correspond to eqs. (1) 
to (6) in the text. The allowed parameter space for PSR 
B0943+10 is represented by the gray area. The two points 
with co-ordinates (T^-g , h) marked by * and • represent 
two modelled cases (1.5, 7.2) and (1.9, 6.7), respectively, 
described in section 3 and presented in Figs. 2 and 3. 



attempt to estimate the speed of circumpferential motion 
of each spark Vd = cE±/Bs, where E± ~ 2/S.V/rp (to 
within an uncertainty factor of less than 4) and /SV = 
{2i: / cP)Bsh? is the potential drop within the gap. Thus, 
the speed of E x B circulation Vd = (47r/P)(/i^/rp). In the 
case of PSR B0943+10, h ^ 500 cm (see Fig.3), Tp - 3300 
cm and Vd ^ 850 cm/s. The time interval after which 
each spark completes full (360°) circulation around the 
polar cap perimeter is P3 ~ 2Trrp/vd=0-5P{rp/h)^ , which 
is about 22P in this case. Given the uncertainty in E±, 
this is consistent with to the observed value P3 = 37P 
( which implies Vd = 27rrp/P3 ~ 510 cm/s, close to 810 
cm/s obtained above given the uncertainty in the theoret- 
ical determination of Ej_). 



2.4- Parameter space 

The surface magnetic field Bg ^ Bq = 4.4 x lO^'^ G repre- 
sents the so-called photon splitting level, above which the 
electron-positron plasma cannot be produced (see Zhang 
& Harding 2000 for review). Hence, yet another constraint 
can be imposed for PSR B0943-hl0 

^" 11. (6) 



Thus, the parameter space for possible values of b and TZq 
is marked by the gray area enclosed by lines (6), (5), (4) 
and (2)0 in Fig. 1. Within this area the ICS-dominated 
NTVG can exist, with 20 equidistant sparks circulating 



^ The line (2) in Fig. 1 corresponds to the case when the 
actual value of the efficiency parameter k = Qrad/qheat = 0.5 is 
adopted. For k — 1 the line (2) should be replaced by the line 
(3). Therefore, the region between lines (2) and (3) corresponds 
to 0.5 < fc < 1, while the region above line (3) corresponds to 



around the perimeter of the polar cap in the time interval 
P3 close to 37 pulsar periods. The CR-NTVG represented 
by line (1) correspond to the case k = 0.5 (for k = 1.0 
the limiting line 6 = 20 x TZq'^^ almost coincides with the 
ordinate axis in Fig. 1; see also section 4). While we are not 
able to exclude the possibility of CR-NTVG formation, we 
find modelling of implied magnetic structure extremely 
difficult as we describe in section 3 and 4. As seen from 
Fig. 1, the surface magnetic field Bg exceeds 2 x 10^^ G 
(6 > 3) and the radius of curvature TZ of surface field 
lines is about (1 - 3) x 10^ cm (1< 7^6< 3). The two filled 
symbols (* and •) represent two model examples of the 
surface magnetic field structure which are discussed in the 
next section. 



3. Modelling the surface magnetic field 

We have argued that the observations of regularly drifting 
subpulses in PSR B0943-I-10 imply that ICS-dominated 
NTVG can form in this pulsar with a strong non-dipolar 
surface magnetic field. Further from Fig 1, we notice 
that the two parameters lie in the range 3<6<11 and 
0.33< pc< 1, where we introduce the normalized curvature 
Pc = l/Ti-e- We now examine by means of numerical mod- 
elling whether such a structure of the surface magnetic 
field is conceivable. We assume that the actual surface 
magnetic field is a superposition of the star centered global 
dipole moment d and the crust anchored dipole moment 
m placed at Rm, whose infiuence results in small scale 
deviations of the surface field from the pure global dipole. 
The resultant surface magnetic field = 'Bd+'Bm, where 
Bdl ~ 2d/r^ and |Bto| ^ 2TO/|r — Rm|^, where r is the ra- 
dius vector. For clarity of graphic presentation we use the 
normalization Rq = |R„j|/P and the strength of rn = |m 
is expressed in units of d = |d|. 

The detailed numerical formalism of modelling the ac- 
tual surface magnetic field wit hin the above model was 
developed by Gil ct al. (2002 ) and we now apply it to 
PSR B0943-f 10. This accompanying paper gives details 
of modelling of the surface magnetic field in pulsars (a 
review concerning the crust origin surface magnetic field 
in pulsars is also included in this paper). To obtain the 
equations of resultant open magnetic field lines (along 
which a high unipolar potential drop can develop) we solve 
the system of differential equations d6/dr = Bg/{rBr) 
and d(j)/dr = Bg{rBrSm9), with the initial condition 
Bmir = 5P) = (since B,„/Bd - im/d){r/\r - Rm\f, 
then for m/d - 10""^ the ratio Bm/Bd ^ lO"'^ for r> 5R). 
Thus, we define a bundle of the open dipolar field lines at 
the altitude r = 5P and then trace the lines of resultant 
magnetic field Bs down to the stellar surface (r = R). 
Consequently, we obtain the boundary of the physical po- 
lar cap as well as the actual surface magnetic field above 



k = 1 (for more detailed discussion of the parameter k see Gil 
(2002), in preparation). 
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Fig. 2. Structure of the surface magnetic field Bs for 
a superposition of the global star centered dipolar mo- 
ment and crust anchored dipole moment described by: 
m = 2.85 X 10"'' d and Rq = 0.967 R (7k-) and m = 
1.73 X 10"'' d and Rq = 0.975 R (•), respectively. The hor- 
izontal axis is labelled by the azimuthal angle (magnetic 
CO- latitude). The vertical axis is labelled by the normal- 
ized surface magnetic field b = Bg/ B^i in the upper panel, 
and by the normalized altitude z — r/Rin the lower panel. 
The angular size of the canonical (^^-dipolar) and the ac- 
tual (0s-non-dipolar) polar cap is marked by dotted and 
dashed vertical lines, respectively. The converging mag- 
netic field lines presented in the lower panel correspond to 
both cases, but, given the proximity of the two structures, 
cannot be distinguished within the graphic resolution. 



scissa is labelled by the azimuthal angle 6 (magnetic co- 
latitude), which measures the polar cap radius in radi- 
ans. For a purely dipolar magnetic field Vp ^ R ■ sinOd ~ 
2Q4p-i/2 which corresponds to \dd\ ^ 0.01 radians 
(marked by two vertical dotted lines). The actual polar 
cap is much narrower, with j^^l ~ 0.004 radians (marked 
by two vertical dashed lines) or ~ OArd- Near the last 
open magnetic field lines (thick lines in the lower panel) 
b = 6.7 (•) and b = 7.2 (★), respectively. 

Fig. 3 shows the average normalized curvature pc — 
1/TZc of the two outermost (thick) magnetic field lines 
shown in Fig. 2. The horizontal axis is labelled by C = 
where z — r/R (Fig. 2). The dashed vertical line corre- 
sponds to the characteristic gap height h w 500 cm. As 
follows from the parameter space marked by the gray area 
in Fig. 1, the normalized radius of curvature 1< 7?.6< 3 and 
thus the normalized curvature 0.33< Pr < 1. Therefore, con- 
ditions for the formation of the ICS-NTVG are satisfied 
within the gray area marked in Fig. 3. This corresponds to 
an annular region at the periphery of the polar cap (with 
size comparable with the gap height hies ^ 0.15rp), and 
can accommodate only one ring of 20 drifting sparks (with 
T) ^ hies)- Within the central region inside this ring, the 
polar gap with sparking discharges cannot develop. Thus, 
either the central part of the beam of PSR B0943-I-10 is 
radio-quiet (hollow-cone) or its radio emission is driven 
by the stationary free-flow (e.g. Sharleman et al. 1978) 
instead of nonstationary spark-associated flow (RS75). 

As it can be seen from Fig. 3, the crust-origin magnetic 
field decays rapidly with altitude and at C > 0.3 (z > 1.3) 
the normalized curvature pc — 0.007 {TZq ~ 150), sug- 
gesting almost for purely dipolar field. Moreover, the sign 
of Pc in the region of dipolar field lines is negative, while 
within the gap (and its vicinity) pc is positive. Thus, the 
geometry of the surface magnetic fields lines within the 
gap is converging, implying that the radius of curvature 
is directed towards the pole (opposite to the case of di- 
verging dipolar field lines). Therefore, the magnetic field 
driven spark discharges will tend to occupy the peripheral 
areas of the polar cap, instead of moving rapidly towards 



the polar cap. It follows from the symmetry suggested by the pole (Cheng fc Rudcrman, 1977; Filippenko fc Rad 



the observed patterns of drifting subpulses and their spec- hakrishnan, 1982). This should have a stabihzing effect 



tral analysis (DR99, DROl) that the local dipole m should 
lie in the fiducial plane containing both the pulsar spin 
axis and the magnetic axis. Without loss of generality, we 
placed for simplicity the local dipole axis at the polar cap 
center. Since Bs > Bd in PSR B0943-I-10 then in order 
to increase Bd to values exceeding 10^^ G the magnetic 
moments d and m must have the same polarity. 

Fig. 2 presents results of calculations for two cases: 
m = 2.1 X 10~^d and Rq — 0.967i? (marked as ★), and 
m = 1.0 X lO-'^d and Rq = 0.975i? (marked as •). The 
lower panel presents the actual surface magnetic fields 
(as a function of normalized altitude z = r / R) for both 
cases, however, the proximity of the structures makes it 
difficult to distinguish them in the plotted scale. The ab- 



on the arrangement of 20 sparks circulating around the 
perimeter of the polar cap in PSR 609434-10. 

4. Summary 

In this paper we explored consequences of the vac- 
uum gap model interpretation of drifting subpulses in 
PSR B0943-I-10, where 20 sparks move circumferentially 
around the perimeter of the polar cap, each complet- 
ing one circulation in 37 pulsar periods (DR99, DROl). 
Both the characteristic spark dimension 2? and mutual 
distance d between adjacent sparks are assumed to be 
about the polar gap height h. The formation of the vac- 
uum gap requires extremely strong non-dipolar surface 



6 



Gil, Melikidze, Mitra.: Vacuum gap model for PSR B0943+10 



0.01 
0.00 
-0.01 





' ' "1 1 ' ' ' ' 


' ' "1 

• = 0.97 ■ 




- 


/ • \ 


m = 0.0001 
ft = 7.18 
91=1.50 

6 






r /;i-7.89 




/ / \ 


- 452 cm 




^ -kR =0.967 






-^"^ m = 0.00021 






h = 6.72 






% =1.94 


W * ' 




rjh = 6.60 






h — 558 cm 






■ ■ ■ ■! 1 ■ ■ ■ ■ 


■ ■ ■ ■! 



10 10 ^ 10 

Fig. 3. Plot of the normalized curvature pc averaged for 
the outermost open magnetic field lines (thick lines in 
Fig. 2) versus the normalized altitude expressed in the 
form = z — 1, where z ~ r/ R (Figs. 2 and 3). The param- 
eters of the modelled points marked by 7k- and • are given. 
The gray area indicates values of pc for which the ICS- 
NTVG can form in PSR B0943+10. The vertical dashed 
line represents the VG height h ^ 500 cm. 



magnetic (GMOl). We considered both the CR and res- 
onant ICS seed photons as sources of electron-positron 
pairs and determined the parameter space for the sur- 
face magnetic field structure in each case. For the CR- 
NTVG {k = 0.5) the surface magnetic field strength 
i?s > 2 X 10^'^ G and the radius of curvature of surface field 
lines 6 x 10'*cm < 7^ < 12 x lO^cm, while for the resonant 
ICS-NTVG Bs>2x lO^^ G and lO^cm < 7^ < 3 x lO^cm 
(in both cases Bg < Bq ^ 4.4 x 10^^ G). Our modelling 
favours ICS-NTVG formation (refer Fig. 1) although CR- 
NTVG cannot be excluded. The ICS photon driven vac- 
uum gap supported by the magnetic field structure de- 
termined by this parameter space guarantees a system of 
20 sparks with V ^ d hies, circulating around the 
perimeter of the polar cap by means of the E x B drift in 
about 37 pulsar periods. This phenomenon is most prob- 
ably observed in the form of drifting subpulses in PSR 
B0943-fl0. 

We modelled the magnetic field structure using a sim- 
ple and quite general model in which the actual surface 
magnetic field at the polar cap is a superposition of global 
core-anchored dipole (inferred from pulsar spin down) 
and local crust-anchored dipole B„i ( Gil ct al., 2002| ). 
Since Ba = 6.4 x 10i2(p.p)i/2G = 4 x 10^^ G, then in or- 
der to obtain B^ ~ (2 3) x 10^'^ G one needs to assume 
B,ri ^ Bd and the same polarity of both components. 
Following the symmetry suggested by observed patterns 
of drifting subpulses in PSR B0943-I-10, we placed the lo- 



cal dipole axis at the polar cap center. Two examples of 
model solutions (6, TV) are shown in Fig. 1 marked as * 
and • within the gray area. More details concerning these 
two solutions are presented in Figs. 2 and 3. Our mod- 
elling yields a number of important conclusions: (i) The 
conditions for the formation of the ICS-NTVG are sat- 
isfied only at peripheral ring-like region of the polar cap 
which can just accommodate a system of 20 E x B drifting 
sparks with V ^ d ^ hjcs- (ii) The surface magnetic field 
lines within the actual gap are converging, which stabilizes 
the E X B drifting sparks by preventing them from rush- 
ing towards the pole (as in the case of diverging dipolar 
field), (iii) No model solutions could be obtained close to 
or above the fine (3) in Fig. 1 (especially for TZq > 1.5), 
corresponding to the efficiency parameter^ k = 1. This 
means that k< 0.8 (• in Fig. 1), at least in PSR B0943-f 10. 
(iv) We find no model solutions in the range 6 ~ 7 to 10 
and 7^6 0.06 to 0.12, corresponding to the CR-NTVG 
parameter space (not shown in Fig. 1). This favors the 
ICS-NTVG in PSR B0943-f 10. We emphasize that conclu- 
sions (ii)-(iv) result from a specific dependence of b on TZi, 
following from central position of the local surface dipole 
Bm, which seems to be justified in this case by apparent 
symmetries inferred from the observational data of PSR 
B0943-M0. 

The sparks inherent to the vacuum gap model not only 
provide a natural explanation for drifting subpulses but 
they also imply a mechanism of their radio emission. Me- 
likidze et al. (2000) proposed a spark-associated model for 
coherent pulsar radio emission by means of curvature ra- 
diation of charged relativistic solitons and we will apply 
this model to the ICS-NTVG in PSR B0943+10 in the 
forthcoming paper. 
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